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Abstract 

A fluctuation law of the energy in freely-decaying, homogeneous and isotropic turbu- 
lence is derived within standard closure hypotheses for 3D incompressible flow. In partic- 
ular, a fluctuation-dissipation relation is derived which relates the strength of a stochastic 
backscatter term in the energy decay equation to the mean of the energy dissipation rate. 
The theory is based on the so-called "effective action" of the energy history and illustrates 
a Rayleigh-Ritz method recently developed to evaluate the effective action approximately 
within probability density-function (PDF) closures. These effective actions generalize the 
Onsager-Machlup action of nonequilibrium statistical mechanics to turbulent flow. They 
yield detailed, concrete predictions for fluctuations, such as multi-time correlation functions 
of arbitrary order, which cannot be obtained by direct PDF methods. They also characterize 
the mean histories by a variational principle. 



1 Introduction 



We consider here the problem of fluctuations of the energy in high Reynolds number turbulence 
decay. The mean energy decay in homogeneous and isotropic turbulence has been the subject 
of many classic investigations. A rather thorough review is contained in [jy]. Section 16. von 
Karman and Howarth Q first derived a power decay law for the mean energy, K*(t) oc (t — to)~ n , 
by means of a hypothesis of complete self-preservation of the spectrum. The fundamental paper 
of Kolmogorov (3| rederived that result, with a precise prediction for the exponent, n = ^y. 
Kolomogorov's original argument assumed, however, the conservation of the Loitsyansky invari- 
ant, which was later called into question by Proudman and Reid Nevertheless, Kolmogorov's 
basic argument may still be carried through under a weaker hypothesis, a "principle of perma- 
nence of the large eddies". This now-standard theory has been discussed in several books and 
reviews: [p|, 0, and Q. According to this theory, the decay exponent n is dependent on the 
initial-data, through the power of the low-wavenumber part of the spectrum. 

Our interest here is in the fluctuations of the energy history during the decay, including 
joint multi-time statistics. The main results have been briefly announced elsewhere ||. Our 
analysis is based on a general approach to fluctuations in irreversible processes, proposed by 
Onsager || and developed in detail by Onsager and Machlup In this method, an "action 
functional" is employed which measures directly the probability of observing a given history as 
a fluctuation event. In particular, the most probable history minimizes this action functional. 
In systems close to thermal equilibrium, there is a standard fluctuation-dissipation relation for 
molecular noise, so that the Onsager-Machlup action has there the physical interpretation of a 
"dissipation function". Onsager's variational principle reduces near equilibrium to a "principle 
of least dissipation" , generalizing the well-known hydrodynamic principle of Rayleigh. 

In its original form, however, Onsager's principle was restricted to weakly noisy systems 
and could not be applied to turbulence, where fluctuations are large. Recently we have pro- 



posed a generalization which applies as well to strongly noisy systems [11, 12 ] . The variational 
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functionals in this theory, or "effective actions" , have experimental consequences for turbulence 
fluctuations and are subject to realizability conditions which arise from positivity of the under- 
lying statistical distributions. For each random variable Z(t) in the flow (where Z may represent 
a velocity at a chosen point, a pressure, a turbulent energy, etc.) there is a corresponding ef- 
fective action T[z], which is a functional of the whole time-history {z(t) : to < t < +00} of the 
variable. The realizability conditions on this action function are (i) that it be nonnegative , 
r[z] > 0, (ii) that it have the ensemble mean z(t) as its unique minimum T[z] = 0, and (Hi) that 
it be convex, Ar[zi] + (1 — A)!?^] > T[\z± + (1 — A)^], < A < 1. As a consequence, the mean 
value z(t) is characterized by a principle of least effective action. Like Onsager's action, this 
functional directly measures the probability of fluctuations of the sample histories away from 
the mean history. The effective action also serves as a generating functional for (irreducible) 
multitime correlation functions of the considered random variable. 

To make the effective action into a practical working tool, efficient and economical approx- 
imation procedures are required to calculate it. In ][ll], 12] we have demonstrated one such 
scheme, a Rayleigh-Ritz variational method inspired by the similar ones already extensively 
used in quantum theory. This variational method is designed to be used in conjunction with 
probability density-function (PDF) closures, such as mapping closures |l3|, |l4|], generalized 



Langevin models |15|, 16], etc. Any reasonable guess of the turbulence statistics may be input 
into the variational method to yield approximations of the effective actions. By this means, pre- 
dictions are obtained for multi-time statistics which are not obtainable by direct PDF methods. 
The additional information about fluctuations has been found to be very useful in evaluating 
the reliability of PDF closures for practical modelling purposes § . 

The contents of this paper are as follows: In Section 2 we very briefly review the standard 
theory of the mean energy decay in high Reynolds number turbulence and, in particular, we 
cast it in the form of a PDF-based moment closure. In Section 3 we evaluate the effective 
action within the standard theoretical hypotheses, by means of the Rayleigh-Ritz algorithm. 
The realizability of the approximate effective action is verified in the small-fluctuation regime 



by means of a Langevin dynamics for the turbulent energy and a fluctuation-dissipation relation 
is derived for the strength of the stochastic noise term. In Section 4 we discuss some of the 
testable consequences of the theory. In particular, the prediction for the 2-time correlation of 
the turbulent energy is given. Also, the direct empirical significance of the effective action is 
discussed, in terms of fluctuations in iV-sample ensemble averages. 



2 Review of Theory for the Mean Decay 

We outline here the standard theory of mean energy decay in a freely-decaying homogeneous and 
isotropic turbulence with random initial data at high Reynolds number, following the reviews 
in [||, ^, [?]]. For convenience, we assume a model energy spectrum 



E(k,t) 



Ak m k < k L (t) 



ae 



/ 3 (t)k' 5 / 3 k L (t) <k< k d (t) (2.1) 



k > k d (t) 

which has been adopted in some previous studies [17, ||]. Such a spectrum may certainly be 



taken at time t = to for the initial velocity statistics. We are assuming as well that there is a 
permanent form of the spectrum, according to which the spectral shape is unchanged in time 
except through its dependence on the parameters e(t), ki(t) and k d (t). Note that the spectrum 
is not self-preserving, or self-similar, in the usual sense discussed in [jlj], which would imply that 
it have the form E(k, t) = ae 2//3 (t)k~ 5 / 3 f \k£(t)) for some length-scale £(t). In fact, the spectrum 
contains two distinct length-scales, the integral or outer scale L(t) = kj^it) and the dissipation 
or inner scale rj(t) = k^ 1 ^). Certain features of the above model are crude caricatures of reality. 
For example, the spectrum should not vanish for k > kd(t) at any time t > to, even if it did 
so initially. However, the spectrum should always show some rapid exponential decay in the 
far dissipation range. It may be easily checked that such refinements do not change any of 
the results below. The important assumption has to do with the low-wavenumber part of the 
spectrum. It was found by Proudman and Reid from the quasinormal closure Ej that there is 



a backscatter term ~ fc 4 in the energy transfer T(k,t). Hence, as long as m < 4 one expects 
that there is a permanence of the low-wavenumber spectrum, in the sense that the power-law 
k m and its coefficient A remain unchanged in time. On the other hand, if m > 4 initially, 
then it is expected that the spectrum with m = 4 will be established at positive times and, 
if m = 4 initially, then the low-wavenumber spectrum will remain of the same form with a 
time-dependent coefficient A(t). Here we always consider m < 4, so that the "permanence of 
large-eddies" should hold. For finiteness of the total energy, m > — 1 must also be imposed 
and, in fact, we usually take m > so that the spectrum decreases asymptotically at very low 
wavenumbers. For grid-generated turbulence it is has been inferred that m ~ 2 [17]. 

The mean decay law can be derived very simply from these hypotheses. One relation ki(t) = 

3 

(yf e2/,3 (^)) 3m+5 * s i m P ose d on the spectral parameters by requiring continuity at k = k^t). An 
additional constraint is obtained at high Reynolds number by evaluating the dissipation rate 

/•oo 

e(t) = 2u / k 2 E(k,t)dk, (2.2) 
Jo 

which, for k L (t) < k d (t), leads to k d (t) = (J^J e 1/4 (t). Only one independent parameter 
is left, which may be taken to be the integral E(t) = Jq°° dkE(k, t) representing mean energy 
at time t. For the above form of the spectrum it is not hard to show that at high Reynolds 
number, when ki(t) <C kd(t), the dissipation is given as 

e(t) = A m EP(t) (2.3) 

3m+5 x 

with A^ 1 = o?l 2 [z^i + I) 2m+2 A~^+i and p = f^p§. Thus, employing the Navier-Stokes 
equation via its energy-balance, one obtains the closed moment equation 

E(t) = -A m EP(t). (2.4) 

This generalizes the decay equation (29) in the paper of Kolmogorov. Its solution with initial 
condition E(to) = Kq gives a prediction for the energy decay law, in the form 

/+ — t*\ ~ n 

K * {t) = Ko (-an (2 - 5) 



with n = 2ffid|. Here At 



A m {p-l)K 



p-i 



-l 

is a constant with units of time, determined 



by the initial mean energy Kq, and £q = to — At is a "virtual time-origin". 

This simple theory may be cast into the form of a PDF closure by assuming as an Ansatz at 
all times t > to a Gaussian random velocity field with zero mean and with spectrum E(k, t) given 
by Eq. Q2.ip above. The assumption of Gaussian statistics was not used in previous works. It is 
not necessary here either, but it makes simpler the analytical labor in applying the Rayleigh- 
Ritz algorithm. In fact, the Gaussian Ansatz will only be used to evaluate averages of 1-point 



velocity moments, which it is known have statistics in actual fact close to Gaussian. See [18 
Ch.VIII. It will be shown later that the use of a non-Gaussian Ansatz would change only a 
constant in the final results. The model spectrum contains one free parameter, which may be 
taken to be the energy per mass E(t) in the Ansatz, or, what is the same, its mean value of 
the quadratic velocity-moment functional iT(r;v) = ^v 2 (r) at some chosen space-point r. By 
statistical homogeneity, the mean value is independent of this choice. The time-dependence 
of E(t) is then determined by projecting the Navier-Stokes dynamics onto this single moment 
function: 

E(t) = (tfK(r)) E(t) (2.6) 

where 

£ = " £ / J^j [(-(v(r)- V) Wi (r) - V lP (r) + v A Wi (r)) • • •] (2.7) 

is the Liouville operator which generates the evolution of PDF's for the Navier-Stokes dynamics, 
C) is the adjoint operator which generates the evolution of observables, and (■ ■ -)e{€) denotes 
average with respect to the model Gaussian velocity with energy E(t). It is easy to see that 
this prescription to determine the time-dependence leads to 

E{t) = -e{t) (2.8) 

which, using Eq.Q2.3|), is clearly equivalent to the one above. However, putting the analysis 



into this form allows us to apply the Rayleigh-Ritz method of [[11], 12] to evaluate the effective 
actions. 
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3 Calculation of the Action 

We calculate here the effective action T[K] of the energy history K(r, t; v) = ^v 2 (r, t). For each 
time t, K(r, t) is a functional on phase space via its dependence on the random initial data v(r) 
at t = to °f the Navier-Stokes solution v(r, t). K(t) is a possible value of this random variable, 
i.e. a numerical time-history. According to the theorem established in [11, the effective 
action is characterized as the stationary point of the "nonequilibrium action functional" 

T[A,p}= / dt(A(t),(d t -C)p(t)) (3-1) 



varied over arbitrary left and right "trial functionals" ^4[v;t] and />[v;i]. (In [12] these were 
denoted ^ L ,^f R , respectively; the hat is used here to denote functionals on the phase-space of 
velocity fields). The variations are performed subject to the constraints of unit overlap 

(A(t),p(t)) = l (3.2) 

and fixed mean (of K(r) not of K(r,t)\) 

(A(t),K(r)p(t))=K(t), (3.3) 

with the initial condition 

p[v;t } = PoH (3-4) 
where Po is the initial Gaussian distribution at t = to, and with the final condition 

i[v;+oo] = l. (3.5) 

Note that {A,p) = J T>vA[v]p[v]. The trial functional p(t) should be taken to vary over 
the space of all probability distributions, while A(t) is varied over the space of all bounded 
observables. In this case the constraints become, more simply, 

(Mt))m = 1 ( 3 - 6 ) 

and 

{A(t)K) m =K(t), (3.7) 



in which (• • -)p(t) denotes the average over the distribution p(t). Likewise, 

T[A, p} = - dt {(d t + C ] )A(t)) m (3.8) 

J to 

is a generally more useful expression for the nonequilibrium action. 

To obtain the exact effective action, trial functionals should be varied over the full spaces. 
However, within the Gaussian Ansatz above, the variation is taken over a restricted class of 
trial functionals. The right functional is just the Gaussian PDF itself: 



P[v;t] = ^yexp 



(3.9) 



with the isotropic spectral tensor 

W) = f^(%-^) (3.10) 
in which E(k,t) is the scalar spectrum of Eq.(2.1) for a variable total energy E(t). M(t) is 



the normalization factor guaranteeing total probability equal to unity. The left trial functional 
within the Gaussian PDF closure is chosen from among arbitrary linear combinations of the 
moment functional K[r; v], which appeared in the closure, and the constant functional = 1: 

i[v; t] = ao(t)l + ai(t)K[r; v]. (3.11) 

The variable functions of time, E(t),ao(t),ai(t) are the trial parameters of the variational 
calculation. 

However, because of the two constraints, Eqs.(^),( |3.7| ), only one of these parameters is 
independent. We shall take it to be E(t). The unit overlap condition Eq.(^) requires that 
ao(t) + cti(t)E(t) = 1, or that 

i[v; t] = 1 + a x {t) (if [r; v] - E{t)) (3.12) 



by eliminating ceo(t). Next a\(t) may be eliminated by using the condition Eq.(3.7). Since 



(A(t)K(r)) m = E(t) + ai (t) ((K 2 (r)) E(t) - E 2 (t)) , (3.13) 



the constraint equation is obtained from an easily calculated average over the Gaussian ensemble 

«j(r))s(t) = l$iy 



parametrized by E(t). Using (v i(r)vj(r)) E u\ = ^5ijE(t), this average is found to be 



(K 2 (r)) m = \E\t). (3.14) 

Note that evaluation of this 1-point moment is the only place where Gaussian statistics is 
employed in the whole calculation. From the imposed condition Eq.([T7]) we then obtain that 

ai {t) = \E-\t)\K{t)-E{t)]. (3.15) 

The action may now be approximated as 



;>oo 

T # [K;E] = - dt ((d t + tf)A(t)) E{t) 
Jt 

= J°° dt ^-^{K -E{t)) E(t) + ai {t)(E{t) + {e) E{t) \ 
= / dt ai (t) E(t)+A m EP(t) 

Jt 1 J 

= -/ dtET 2 (t)[K(t)-E(t)] E(t) + A m EP(t) , (3.16) 

2 Jta 1 J 



in which E(t) remains as the only trial parameter. We wrote as e(r) = | J2ij(di v j( r ) ~\~dj v i{ r )) 2 
the local energy dissipation rate and noted its average from Eq. (|2.3|) as (i) E M = ^■mE p (t). By 
requiring stationarity of the action under variations of E{t), or, 5T*[K; E]/5E(t) = 0, with 
fixed K(t), it is straightforward to derive the variational equation 

A m EP(t) + K(t) = (p - 2)A m (K(t) - E{t))EV-\t). (3.17) 

For any rational value of p = j, k, I integers, this is a polynomial of degree k in X = E 1 / 1 : 
(p - l)A m X k - (p - 2)A m KX k - 1 + K = 0. For a physically allowable energy history, K{t) > 
and K(t) < 0. Furthermore, p > 1 whenever m > —3. Thus — independent of the sign of 
{p — 2) — the polynomial has one change of sign in its coefficients for any permissable energy 
history. Therefore, it follows from Descartes' rule of signs that there is exactly one positive root 
E{t) for each physical choice of K(t), when p is rational. Because these are dense in the real 



p > 0, E is uniquely defined for all permissable K. Substituting that value into the Eq. (|3.16|) 



above, we obtain the final form of the approximate effective action 

3 ,00 (if (i) + A m EP(t)) (E(t) + A m EP(t)) 

= w^J t , * 1 " (3 - 18> 



in which the -E-dependence is eliminated by inserting the root of Eq. fl3.17 ) as described. 



It is easy to check that, if the approximate action is evaluated at the predicted closure mean 
energy if*(i), then I 1 * [if*] = 0. In fact, using if*(i) = — A m if*(i), the corresponding E*(t) is 
determined from 

A m (E*(t) - Kl{t)) = (p - 2)A m (if*(i) - E.(t)) E^~\t). (3.19) 

The solution of this equation is 

E*(t)=K*(t) & E*(t) = -A m EP(t). (3.20) 

Obviously, substituting these values makes the approximate action vanish identically. It can, in 
fact, be shown that the mean value for any closure is a zero of the approximate action evaluated 



by the Rayleigh-Ritz method within that same closure [11, 12 1. It may even be shown further 
that the mean value is always a stationary point of the action, 5T*[K*]/6K (t) = 0. However, it 
need not be a minimum point, as required by the realizability conditions on the effective action. 

To examine the issue of realizability here, we consider small perturbations K(t) = K*(t) + 
SK(t) from the predicted mean. Because the calculation is straightforward but somewhat 
tedious, we give the details in Appendix I. The final result is that 

3 ,00 (sk{t) + A mP Kr\t)SK{t)) 2 , . 

r - 1*1 - L * 1 SFcij + ° ^ ■ (3 ' 21> 

Note again that the coefficient (p — 1) in front of the action is > as long as m > —3. 
In fact, m > — 1 is required to give a finite energy. Thus, for all permissable values of m, 
the approximate action I 1 * [if] satisfies realizability, at least in a small neighborhood of the 
mean energy history if*(t). One should be cautioned that satisfaction of realizability is only a 
consistency check and cannot guarantee correctness of predictions. Indeed, the same calculation 
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as we made above would carry through exactly for the ID Burgers equation, since the only 
property of the nonlinear dynamics that was used was energy conservation. However, not all of 
the previous results are true for Burgers turbulence. In that case the energy spectrum Eq.( |2.1[ ) 
is not even the correct quasi-equilibrium form but, instead, a k~ 2 spectrum will develop p9| . 
This graphically illustrates that realizability is perfectly compatible with falsity. On the other 
hand, we expect that the approximation to the effective action with the Kolmogorov spectrum 



Eq.(2.1) is qualitatively correct for 3D Navier-Stokes turbulence. Unfortunately, we have not 



so far been able to show that the full action, Eq.(3.18), satisfies all realizability constraints for 
arbitrarily large deviations 5K. 



It may be observed from Eq,( |3.21 ) that the quadratic part of the approximate action has 



precisely the form of an Onsager-Machlup action [10]. Hence, it follows from the work of 
Onsager and Machlup that the same law of fluctuations is realized with the Langevin equation 

SK+(t) + k mP Kr l {t)8K + {t) = (2R,(t)) 1 / 2 r,(t) (3.22) 

obtained by linearization of the energy-decay equation around its solution K* (t) and by addition 
of a white-noise random force rj(t), {t](t) r j{^')) = ~ 0> w hh a coefficient 

R*(t) = 2 ^ P ~ l \ s)K*{t)- (3-23) 

This alternative stochastic representation is equivalent in the sense that all finite distributions 
of 5K + (t) in the above Langevin model agree with those predicted for 5K(t) by the quadratic 
action. We emphasize that this linear Langevin representation is only adequate for the 
smaller fluctuations about the mean and will not be sufficient to describe the larger fluctuations. 
In fact, the quadratic part of the action is only a valid approximation for small deviations 
5K. The predicted decay of the smaller energy fluctuations according to a linearized law 
is similar to the Onsager regression hypothesis for equilibrium fluctuations |J. Likewise, the 



1 Recall from JO], [T^] that r[K] is a generating functional for irreducible multitime correlation functions of 
the energy K(t). See also Section 3. 
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expression Eq.( |3.23| ) is a fluctuation- dissipation relation (FDR) analogous to that in equilibrium. 
The white-noise term on the righthand side of Eq.( |3,22| ) represents a stochastic backscatter 
contribution to the energy evolution and the Eq.( |3.23 ) relates its magnitude to the mean energy 



dissipation rate e*(t). These are testable predictions of the closure hypotheses. We expect that 
the prefactor C = |(p — 1) in the FDR, whose precise value results from the Gaussian Ansatz, 
is correct at least on order of magnitude. For m = 2 its value is C = | ~ 0.56. 

The Gaussian Ansatz is obviously inadequate in one respect, because it fails to capture the 
important non-Gaussian effect of scale energy transfer. Let tl(k) be the usual spectral flux as an 
instantaneous variable in individual realizations, written in terms of a triple product of velocity 
Fourier coefficients. (For example, see [^] for an explicit expression). Then, one expects for 
freely-decaying turbulence in the quasi-steady regime that (JI(k))E(t) = for all inertial- 
range wavenumbers fc^(f) <C k <C kd(t), breaking time-reversal symmetry. However, within 
the Gaussian Ansatz (Il(k))E(t) = 0- This pathology of the Gaussian Ansatz shows up if one 
calculates r*[II], the Gaussian approximation to the effective action of the flux TI(k). In fact, 
r*[n] is the Legendre transform of an approximate cumulant-generating functional A* [if], using 



the notations of [12]. That is, r*[II] = max# (HU — A* [if]). Because the 5th-order moment 
(K (r)tl(k)) ^(t) = in the Gaussian Ansatz, as well as (TL(k))E(t) = 0) h follows by the methods 
discussed in |l^] that A* [if] = for all H. Therefore, the Legendre transform is 

[ o n = o 

r*[n] = <^ (3.24) 
[ +oo n / o 

This result just implies that, within the Gaussian Ansatz, the flux function tl(k) is identically 
zero in every realization and no fluctuations from that value may occur. This is clearly an 
unphysical result of the closure. However, the failure of the Gaussian Ansatz to describe the 
energy transfer is hoped not to drastically affect the result for the effective action T[K] of the 
energy, because the model spectrum Eq.([2.1|) has built-in the correct overall decay rate. 

Some insight into this may be obtained by considering the exact equation for the 2-time 
correlation of the energy fluctuation. The energy density fluctuation 8K(r, t) = K(r, t) — (K(t)) 
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in each individual realization obeys the equation 

SK = -V -[{K + p)w - vV K] -Si. (3.25) 

Here Si = i — (i) is the energy dissipation fluctuation and p is the kinematic pressure. It is 
this equation which is being statistically modeled by the Langevin equation, Eq.(p2|). If this 



model is to be valid for second order statistics, then it must be true that the exact equation 

(sk(t)sk(t )) = -{v-[(k + p) v - uvk}(t)sk(t )) - (Si(t)sk(t )) (3.26) 

coincides with the one derived from the Langevin equation. For t > to this is just 

(Sk + (t)SK + (t )} = ([-L,(t)SK + (t) + (2R,(t)) 1 / 2 7 ] (t)]SK + (t )) 

= -L*(t){SK + (t)SK + (t )). (3.27) 

We have introduced L*(f) = h m pK*~ l (t) and noted that the white-noise force is uncorrelated 
with earlier values of the energy fluctuation. In order to coincide with this model equation, it 
is clear that the first term due to space transport on the LHS of Eq.( 3.26| ) should be negligible, 



i.e. (V-[(k + p)v — u'V k](t)Sk(to)} 0. This is plausible, because the space transport term 
is rapidly varying in time and thus decorrelated with the energy fluctuation at earlier time. 
It is for this reason that such higher moments (of 4th and 5th order in velocity) were never 
explicitly modeled in our analysis, although they are implicitly represented by the white-noise 
term in the Langevin equation. The remaining term in the exact equation coincides with that 
in the model equation, if it is further assumed that 

(Si(t)\k(s),s < t) m -L*(t)5K(t). (3.28) 

That is, the conditional expectation of the dissipation fluctuation given the value of the energy 
density over the entire past should be obtained by linearizing the expression for mean dissipation 
in the closure and evaluating it at the present value of the energy fluctuation in the given 
realization. This conditional relation is assumed to hold when the energies k(t) are small 
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deviations from the mean value {K(t)) ~ K*(t). This formula also has considerable plausibility: 
it is a formal statement of the "regression hypothesis" on small fluctuations. If this relation is 
used to eliminate 5i(t) in Eq.( |3.26D , then an equation of the same form as Eq.( |3.27 ) is obtained: 



(5K(t)5K(tQ)) ~ —L*(t)(5K(t)5K(to)). This should make more transparent the nature of the 
approximations in the Rayleigh-Ritz calculation at the level of closure considered here. 

In fact, it is possible by such arguments to completely "rederive" the Langevin model. If 
one accepts (i) the hypothesis that the rapidly changing terms are correctly modeled by a white 
noise, i.e. 

- V-[(K + p)v-vVK](t) - [5£(t) - (8e{t)\K{s),s < t)} « ^(i)) 1 / 2 /?(*), (3.29) 



and (ii) the "regression hypothesis" in Eq.( 3.28 ), then the exact equation Eq.( 3.25| ) reduces to 



the Langevin model Eq. (|3.22j ). A Kolmogorov-style dimensional analysis would yield for the 
noise strength = Ce*(t)i^*(t) with C some universal constant, to be determined. The 

value of this constant C = |(p— 1) resulting from the variational calculation with the Gaussian 
Ansatz will furthermore be shown below to be the unique choice to recover the relation Eq.( [4.5[ ). 
Since this relation is exact when single-point statistics of velocity are Gaussian — which is known 
to be a very good approximation — the Langevin model can be entirely motivated by intuitive 
considerations. The variational calculation is a systematic analytical procedure yielding the 
same Langevin model, but also only approximate. The two derivations are therefore very 
complementary. 

Improvements of the Gaussian Ansatz are likely to lead to better results for the effective 



actions. For example, the "synthetic turbulence" models of [21| are random velocity fields which 
contain the correct energy transfer and also some of the intermittency effects of real turbulent 
velocities. Using such statistical models within our Rayleigh-Ritz scheme should lead not only 
to a qualitatively correct result for r[II] but also to quantitatively better results for T[if]. In 
such improved closures new "test functionals" in addition to the quadratic moment-functional 
K(r; v) must be considered to determine the time-dependence of the additional free parameters 
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in the statistical Ansatz. For example, the energy flux variable (a triple moment-functional) 
would be a natural variable to add to the closure. The choice of the "test functionals" is an 
equally important element of the closure as is the choice of the model velocity statistics. We 
emphasize again that our results above for T[K] depend very little on the choice of the Gaussian 
statistics. For any model statistics with the mean energy E(t) as the only free parameter and 
with K(r; v) the corresponding "test functional", results very similar to those above will follow. 
In that general setting a result (K 2 (r)) E ^ = BE 2 (t) will hold by dimensional analysis, for some 
constant B, replacing Eq,( [3.14] ). This means that the formula Eq.( 3.1§| ) for the approximate 
action will still hold, with the factor | in front simply replaced by another number D = 1/(1?— 1) 
of order unity. Only the value B = | depends upon the Gaussian Ansatz. By employing 
improved closures one may hope to derive from first principles such theoretical features as the 
"permanence of large eddies" for m < 4. Because the Rayleigh-Ritz algorithm is a convergent 
approximation scheme for the true effective actions, systematic improvement of the closures 
will lead to a refined description of the turbulent dynamics. 

4 Testing the Theory 

The previous theory has testable consequences for turbulent energy fluctuations. The most 
likely experimental situation for such checks is grid turbulence, which well approximates a 
homogeneous, isotropic, decaying turbulence. In principle, it would be possible to make an 
experiment by measuring the velocity at a single point r in the frame of mean downstream 
motion. Constructing from this the energy history K{t) = iv 2 (r, i) and compiling an ensemble 



of realizations, one may, as in [17, 22], compute various multi-time statistics to compare with 
the theory. 

The most familiar such statistics are the r-time correlation functions (K(t\) ■ ■ ■ K(t r )). 
These may be derived directly from the effective action r[_RT]. In fact, by taking r functional 
derivatives of the action, evaluated at the mean value, the irreducible r-time correlators are 
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obtained: 



{K{t 1 )---K(t r )) i 



5 r T[K] 



8K{t x )---5K{t r 



(4.1) 



K=K, 



For this result, for the definition of irreducible correlators and their relation to the connected 



correlators (or cumulants), see any text in quantum field theory, e.g. |p3f |. Section 6.2.2 or [24], 
Section 10.2. We only note here that the irreducible 2-time correlator, or (K(ti)K(t2)Y" , is 
the inverse operator kernel of the connected 2-time function (covariance) {K(ti)K(t2)) con = 
(SK^SKih)}, i.e. 

' ds {K(t)K(s)Y rr (5K(s)5K(t')) = 5(t - t'). (4.2 



Relations between higher-order irreducible and connected correlators are obtained by taking 



further functional derivatives of this relation with respect to K. See p3| , 24] 



It is very easy to obtain the variational approximation (K(t)K(s))" T from Eq.(4.1) and the 
quadratic part of the Gaussian effective action, Eq.( 3.2l| ). Taking its inverse, the covariance 
(5K(t)8K(t'))* is then evaluated as 



(6K(t)5K(t'))* = exp 



t ft' 
ds L*(s) - ds L*(s) 

to J to 



rmm{t,t'} 

+2 / ds R* (s) exp 

Jto 



(SK ) 



dr LJr) 



dr L*{r) 



(4.3) 



in our theory. We have again written L*(t) = A m pK* x (t). These calculations are outlined 
in Appendix II. In the same way, by taking an arbitrary number r of functional derivatives in 



Eq.(4.1), all correlations of any finite order are obtainable from the effective action. However, 
we do not pursue the general calculation here. 

To cast the theoretical results into a form that may be compared with experiment, we insert 
the mean decay law K*(t) from Eq.([2.5|) into Eq.( [4.3[) and perform the integrals. For the 2-time 
covariance of the turbulent energy this calculation is straightforward and the prediction is: 



(6K(t)6K(t'))i 



t _ t *y(n+l) /f/_t*N-("+l) 



At 



(«5if f + 3^ 2 



At 

t 



— t* 

mm •'O 

~~ At 



(4.4) 
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with i m i n = min{i, t'}. The notations are the same as for the mean decay law. It should be noted 
that the first term oc (<5-K"o) 2 corresponds to decay of an initial energy fluctuation 5Kq. The 
second term oc Kq represents the new fluctuations generated by the internal turbulence noise, 
through the stochastic backscatter dynamics. As a consequence of that term, the long-time rms 
value of the energy, K rms (t) = ((SK(t)) )* , evolves to a constant level with respect to the 
mean energy K*(t): 

lim «'*('»'>• _ I (4.5) 

The limiting value of | is what would occur for an asymptotic Gaussian statistics of the 1- 
point velocity variable. If we had adopted a non-Gaussian Ansatz in our calculation, then the 
predicted limiting value would have been 1/D = (B — 1). Since any value of the constant can 
be accommodated by an appropriate such Ansatz, it is not so important to the theory which 
particular constant is correct (although the Gaussian value is expected to be quite accurate). 
What would falsify the present theory would be an experimental finding that the functional 
form in Eq.(4.4) was wrong, for any possible choice of the constant 1/D replacing |. We should 
emphasize, however, that the standard theory for the mean energy decay law K*(t) in Eq.( |2.5[ ) 
is in agreement with present experiments, with a value of n near 1.2. In all cases studied so far, 
PDF Ansdtze adequate to predict the mean values of selected variables have also, employed in 
our Rayleigh-Ritz method, yielded good predictions for the fluctuations of those variables near 
the means. See M. We therefore expect the prediction in Eq.(fOI) to be reasonably accurate. 



The previous result for r = 2, the covariance in Eq.(4.3), may be obtained as well from the 
Langevin equation, Eq.( |3.22| ). However, it must be emphasized that only the small fluctuations 
of the energy, with 5K(t) <C K*(t), are expected to be distributed according to that linearized 
equation. Because correlation functions will get sizable contributions from the larger fluctua- 
tions, for which the linear law breaks down, it would not be appropriate to compare general 
r-time correlation functions of K{t) obtained from the linear theory with experiment. It only 
happens for r = 2 that the linear Langevin equation and the full (nonlinear) effective action 
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yield the same predictions. 

It is possible to give both the effective action functional and the linear Langevin equation 
a direct empirical significance in grid turbulence. This is based upon the standard device of 
making N independent trials to calculate the averages from experiment. Indeed, performing 
the same decay experiment N times independently, one usually considers an empirical mean 
history 

— 1 JL » 

KN(t) = -Y, K ^) ( 4 - 6 ) 
i=i 

formed from the realizations Ki(t), i = 1,2, ...,N, of the N different samples, i The effective 
action T[K] measures the probability for the empirical mean Kjf(t) (which is a random quantity 
at finite N) to take on a value very different from the true ensemble-average (K(t)). More 
precisely, 

Prob ({K N {t) « K(t) :t <t< +oo}) ~ exp (-NT[K]) . (4.7) 

Thus, the probability to observe Kpfit) taking any value K (t) other than the true ensemble- 
average (K(t)) is exponentially small in the number of samples N. This is a consequence 
of the famous Cramer theorem on large-deviations of sums of independent random variables 
(e.g. see @, Section 8.6.4 and references therein). Put another way, the dimensionless quantity 
l/rfif] gives an estimate of the number N of additional independent samples required to reduce 
by e-fold the probability of the fluctuation value K in the empirical average K^. In principle, 
therefore, the effective action T[K] is itself directly measurable in grid turbulence, by assembling 
a histogram of observed histories K n (t) and determining the decay rate of the probabilities for 
large N. However, this would not be feasible with a reasonable number of independent samples 
N except for histories K(t) sufficiently near the mean history. 

The linear Langevin model Eq.( 3.22j ) is more restricted in its validity, since, as has been 



stressed, it is equivalent to the quadratic action and is adequate only to predict 2nd-order 



2 Assuming that the turbulence is indeed statistically homogeneous, these iV measurements might even be 
taken from points Ti, i = 1, N in the same flow but at separations greater than L, the integral scale, to assure 
statistical independence. 
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statistics of K{t). In particular, a linear Langevin equation can produce only Gaussian multi- 
time statistics for 8K + {i). (Of course, this has nothing to do with the use of a Gaussian PDF 
Ansatz for the velocity field and will be true even if a non-Gaussian Ansatz is employed: see 
Appendix I.) The true statistics of 5K(t) will not be Gaussian at all, e.g. they will be "chi- 
square" if the 1-point velocity itself is Gaussian. Nevertheless, the Langevin equation can also 
be given a direct empirical significance in terms of the independent iV-sample ensemble, based 
upon the central limit theorem. It accurately describes the statistics of the normalized sum 
variable: 

1 * . 

5K N (t) = — - J2(Ki(t) ~ <*(«)»■ (4-8) 



N 1 / 2 . , 
1=1 



This quantity has zero mean and the same covariance as 5Ki(t) = Ki{t) — (Ki(t)) for each inde- 
pendent sample i, i.e. (5K]y(t)SKiy(t')) = (5Ki(t)5Ki(t')) for all i. However, it is furthermore 
a Gaussian variable in the limit of large N, in agreement with the linear Langevin equation. 
Thus, at a large but finite N it is legitimate to compare predictions of the correlations of 5K + (t) 
using the linear Langevin dynamics, Eq.( [3.22| ), with those from experiment for 5KK[{t) at large 
N. Notice that the Gaussian statistics for 5K^{t) in fact result by substituting into T[K] in 
Eq.(l4~7p the value K = (K) + N'^SK. In that case, by expanding in 5K, one obtains 

Prob ({SK N (t) w SK(t) :t <t< +cxo}) ~ exp (-r (2) [5A']) , (4.9) 

where r^f&A] is the quadratic approximation to the exact effective action and terms in the 
exponent of order iV -1 / 2 have been neglected in the large N limit. This is just one of the 
standard proofs of the central limit theorem. The important point here is that it naturally 
accounts why the linear Langevin equation may be appropriate to calculate 2nd-order statistics 
but certainly not higher order. For the latter purpose the full nonlinear action, Eq.( |3.18 ), must 
be used, not just the quadratic part. 
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5 Conclusions 



The main results of this work are as follows: 



(1) We have derived an action functional, Eq.( 3.18| ), for energy histories in decaying ho- 
mogeneous and isotropic turbulence at high Reynolds number, by means of a Rayleigh-Ritz 
calculation using standard closure assumptions. This action generalizes the Onsager-Machlup 
action to fully-developed turbulent flow and characterizes the mean energy history by a varia- 
tional principle. 

(2) We have shown that the quadratic part of the action, Eq.( |3.2l| ), valid for a region 
of small fluctuations sufficiently near the mean, satisfies all required realizability constraints. 
In fact, it is of the Onsager-Machlup form and thus has a stochastic realization by a linear 
Langevin equation for the energy history. The deterministic part of the equation is obtained 
by linearization of the mean decay law and the random part has its strength determined by a 
"fluctuation-dissipation relation" in terms of the mean energy dissipation. 

(3) Testable consequence of the theory are r-time correlation functions of the energy, which 
may be obtained by functional differentiation of the effective action. These are not obtainable 
from the starting hypotheses by direct PDF methods. As an example, the 2-time cumulant, or 
covariance, of the energy history is derived in detail. 

(4) A direct empirical significance of the effective action is given in terms of fluctuation 
probabilities for ensemble averages over N independent samples or ensemble-points. This inter- 
pretation permits the effective action itself to be measured experimentally, at least for arguments 
in path-space sufficiently close to the mean history. 
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6 Appendices 



Appendix I: Quadratic Approximation to the Action 

Let us consider in Eq.( 3.1S| ) small fluctuations 5K,SE: 



K(t) = K*(t) + 5K(t) & E(t) = K*(t) + SE(t), 



(6.1) 



using the fact that E*(t) = K*(t). It is then easy to see from Eq.( |3.18 ) that, up to quadratic 
order, 

-1/ 



oo (SK(t) + A mP Kr\t)5E(t)) (SE(t) + K mP Kr l {t)5E{t) 
dt 



2(p-2)A m 4 ' Kl + \t) 
By a straightforward linearization of the determining Eq. d3.17 ), it follows that 



- 2(p - l)K m KP~\5E - 5K) = 5K + A^K^SK. 
Let us introduce a shorthand notation for the lefthand side of this equation: 

A = 5K + K mP KP- l 5K. 
Using the Eq.( |6.3| ), it is not hard to show that 

6K(t) + K mP Kr\t)8E{t) = „ (P ~ 2 \ A 



and 



SE(t) + k m pKP- l {t)5E{t) = Ia-K 



2(p-l) 



d 



2 "*df \2(p- l)K m KlJ ' 
Substituting these into Eq.( |6.2D above, it follows that, to quadratic order, 



3 


2(p 


-2)A m 




3 


8(p 


-l)Am 




3 


8(p 


- l)A m 



- (p-2) A 



eft 



2(P - 1) *t +1 
A 2 (t) 



A 



2 *dt V2(p- l)A m K p 

3 [°° dt d T ' * ^ 2 

3 l^ +1 (t) 16(p-l) 2 A^7i ^ 

°°eft^. 



A 



(6.2) 



(6.3) 



(6.4) 



(6.5) 



(6.6) 



(6.7) 
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To obtain the last line we used the boundary conditions A(to) = A(oo) = 0, which are required 
by the Eqs.fl3~ip, ( gg) , ( g3| ). It should be obvious that the last line of Eq.fl677D is the same as 
the result, Eq.( |3.2i ), claimed in the text. 

Although the present calculation employed the Gaussian Ansatz, it should be stressed that 
a similar result will hold for more realistic statistical models of the velocity field. In particular, 
the quadratic form of the action does not depend upon the Gaussian assumption, but is simply 
a consequence of the fact that the mean history K*(t) is required to be an absolute minimum. 
Hence, this will be true for any closure model leading to an effection action satisfying the 
realizability conditions. In that case, an expansion in small deviations 5K{t) around the mean 
must necessarily lead to a quadratic expression involving the linearized evolution expression, 
A(i) = 5K+A m pK*~ 1 5K. What will be different for other closures is the coefficient multiplying 
A 2 (t), which correspond to different predictions of the fluctuations around the mean. 
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Appendix II: The Predicted 2-Time Cumulant 

As observed in the text, the irreducible r-time correlations of K(t) can be obtained from func- 
tional derivatives of F[K] at K = K*: see Eq.( [4.1[ ), Equivalently, these irreducible correlators 
can be read off from the Taylor series: 



00 1 r r 

r t X ] = E- dt r T r (t 1 ,---,t r )SK(t 1 )---5K(t r ), 

r=2 T - J J 



(6i 



as (K(ti) ■ ■ ■ K(t r )) lvv = T r (ti, ■ ■ ■ ,t r ). It is thus easy to obtain in the Gaussian approximation 
(K(ti)K(t 2 ))^ T from the quadratic part of the Gaussian action derived in the previous Appendix 
I, written as 

(SKitx) + Ljt l )SK(t 1 )) (6k(t 2 ) + L*(t 2 )6K(t 2 )] 



T*[K} = ^[ dh J dt 2 



5(h-t 2 ). (6.9) 



Thence, 



(^(ii)K(t 2 ))i rr = [-d tl + Ljt x )\ {2R,{t l )y l [8 tl + L.(ti)]5(ti - t 2 ) 
and, from Eq, (|4.2| ), the 2-time cumulant must satisfy 



d tl + L m (ti)] (2R,(t 1 )y 1 [d tl + L*(ti)](6K(ti)6K(t2)) m = 5(h - t 2 ). 



(6.10) 



(6.11) 



To solve this equation, we use the Greens functions 

[-8 + (t,0 = exp 
[d + L^itj') = exp 



ds LJs) 



ds LJs) 



9{t' - t) 
0(t-t'). 



(6.12) 
(6.13) 



which are anti-causal and causal, resp. Performing one integration we obtain, for s <t 2 , 

[d s +Lj^s)}(6k(s)6k(t 2 ))* = 2R^s)[-d s + L^s)]~ 1 5(s-t 2 ) 

= 2i?*(s) exp - / 2 dr L m (r) 9(t 2 - s) 



= G(s;t 2 ). 



A second integration forward from to gives, for t% <t 2 , 



(6.14) 



(SK^dkfa))* = exp 



to 



ds Lj^s) 



f ti 

F{t 2 ) + / ds exp 
Jt 



ti 



dr Ljr) 



G(s;t 2 ). (6.15) 
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The first term on the right is an arbitrary solution of the homogeneous equation. The function 
F(t 2 ) is determined by symmetry and initial conditions as 

rt 2 



F(t 2 ) = exp 



ds L*(s) 



{(5K(t )) ). 



(6.16) 



When this is substituted into Eq.( |6.15 ) along with the definition of G(s; t 2 ) from Eq. Q6.14 ), the 
claimed result, Eq,( |4.3| ) in the text, is obtained. Note that ((6K(to)j ) defines the quantity 
{6K ) 2 in Eq.Q of the text. 

It is easy to check that this result coincides also with (5K + (t)5K + (to)) calculated from the 
Langevin model, Eq,( |3.22 ). However, in general, for correlations of order r > 2, the Langevin 
model is not adequate and the above method of calculation from the effective action must be 
employed. 
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